Cosmological neutrino mass detection: The best probe of neutrino lifetime 
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Future cosmological data may be sensitive to the effects of a finite sum of neutrino masses even 
as small as ~0.06 eV, the lower limit guaranteed by neutrino oscillation experiments. We show that 
a cosmological detection of neutrino mass at that level would improve by many orders of magnitude 
the existing limits on neutrino lifetime, and as a consequence on neutrino secret interactions with 
(quasi-)massless particles as in majoron models. On the other hand, neutrino decay may provide a 
way-out to explain a discrepancy < 0.1 eV between cosmic neutrino bounds and Lab data. 
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PACS numbers: 13.35.Hb, 98.80.-k, 95.30.Cq 



Introduction Recent years have seen an impressive 
improvement on the cosmological constraints to the sum 
of neutrino masses £ = ^2m v (for reviews see [H 
with current limits typically ranging below 1 eV, and 
the most aggressive bounds (but also the most fragile 
ones with respect to unaccounted systematics) already 
at £ < 0.2 eV, 95% C.L. Several forecast analyses 

suggest that cosmological probes will reach in the fu- 
ture an incredible sensitivity to the effects of even a tiny 
mass of the cosmic background neutrinos. In particular, 
cosmic microwave background (CMB) lensing extraction 
may be sensitive to £ ~ 0.035 eV [5|; CMB plus weak 
galaxy lensing with tomography may also push the sen- 
sitivity to £ below the level of ~ 0.05 eV pi, 0], with an 
error as low as ~ 0.013 eV @. Also galaxy cluster sur- 
veys may probe E down to ~ 0.03 eV [8|, and a sensitivity 
down to 0.05 ± 0.015 eV may be reached combining CMB 
with the data from the Square Kilometre Array survey 
of large scale structures [9j. These forecasts show that 
cosmology has a potential sensitivity to neutrino masses 
well below the 0.1 eV level. Of course, the ultimate level 
of the systematics to beat has yet to be reliably estab- 
lished. On the other hand, the synergy between different 
strategies and probes may help to identify the systemat- 
ics, and also to improve over the above-mentioned figures 
of merit. 

The interest of these expectations relies on the fact 
that neutrino oscillation data imply E > 0.06 eV, where 
the minimum E ~ 0.061 ± 0.004 eV is attained for a nor- 
mal hierarchy (NH; values quoted at 2 cr, see [3]). For the 
case of an inverted mass hierarchy (IH), the oscillation 
data imply £ ~ 0.1 eV. In the following, we shall proceed 
under the assumption that cosmological observations will 
be able to probe these Lab predictions. To be defined, we 
shall assume that neutrinos have a hierarchical spectrum 
of either inverted or normal sign, as favored by many the- 
oretical models, including the simplest seesaw ones. A 
degenerate mass scenario is phenomenologically allowed, 
with the existing constraints given by E < 6 — 7 eV if only 
laboratory bounds from tritium endpoint flo| are used, 
or E < 0.2 — 2.0 eV from existing cosmological observa- 
tions, where the range depends on the datasets and priors 
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assumed [l|, The degenerate scenario can be tested 
to some extent independently from cosmological obser- 
vations via future tritium endpoint spectrum ll| and 
(if neutrinos are majorana particles) neutrinoless double 
beta decay [l2T ] experiments. We want to remark, how- 
ever, that our considerations would apply qualitatively 
to a mildly degenerate mass pattern, too. 

The main point of this paper is to motivate that, if a 
positive cosmological mass detection is achieved as ex- 
pected, one will be able to put a remarkably strong con- 
straint on the neutrino lifetime. Note that previous at- 
tention has been paid to the cosmological signatures of 
decaying neutrinos Yet the mass range explored in 
those papers is very large compared to present bounds, 
and the main signature considered was the impact on the 
integrated Sachs- Wolfe effect on the CMB. In our consid- 
erations, the bound comes from the impact that massive 
neutrinos have in the background evolution of the uni- 
verse, in a range of masses where they are relativistic 
well after the CMB decoupling. 

Bounds on neutrino lifetimes are usually quoted in 
terms of the rest-frame lifetime to mass ratio r/m. Given 
a measurement in the time interval t using neutrinos with 
Lab energy E, the naive bound which one can put is 
r/m > t/E. Using then the longest timescale available, 
the universe lifetime to ~ Hq 1 (where Hq is the Hubble 
constant), and the lowest energy neutrinos, the ones of 
the cosmic background which are at least partially non- 
rclativistic, a bound of the order of (m.50 = to/50 meV) 
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is the strongest constraint attainable in principle. This is 
to be compared with the strongest direct bound available 
at present given by the observation of solar MeV neu- 
trinos, of the order of ~ 10~ 4 s/eV [141 ] . and the much 
stronger (but model dependent) bound r/m > 10 11 s/eV, 
which might derive from the observations of diffuse super- 
nova neutrino background [l5l ]. Recently, a bound com- 



parable to those projections, t/to > 4 x 10 



s/eV, 



has been claimed to follow already from the requirement 
that the neutrinos are free-streaming at the time of the 
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photon decoupling, as deduced by precise measurements 
of the CMB acoustic peaks [l6j |. Yet, the robustness of 



this conclusion has been questioned in [17|. We shall see 
that the proposed bound based on cosmological neutrino 
mass detection would be much closer to the maximal the- 
oretical bound of Eq. |T]), thus superseding by several or- 
ders of magnitude the previous ones. More importantly, 
it is not based on a model for secret neutrino interac- 
tions, but on the "observation" of neutrino survival, and 
it applies whatever the final state light particles are. 

The bound In order to estimate the bound on the 
neutrino lifetime from the cosmological observation of the 
neutrino mass, let us recall first how massive neutrinos 
affect cosmological observables. We shall base the follow- 
ing discussion mostly on the treatment given in [jO] , which 
we address for details and further references. The main 
effect is due to the direct or indirect impact of massive 
neutrinos in the background evolution (Friedmann law) 
of the universe. In particular, in the matter epoch (but 
the extension to a late dark energy dominated phase is 
straightforward [H) and assuming stable neutrinos the 
Friedmann equation writes 
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where a dot represents a derivative with respect to the 
conformal time r], a is the scale factor of the universe, 
Gjv is Newton's constant, p m is the average cold dark 
matter (CDM) plus baryon density, and p v is the neu- 
trino energy density; we have defined /„ = p v j (p m + p u ) 
and the second equality in Eq. @ holds to first order 
in /„. Keeping e.g. p m + p u constant, a non-vanishing 
/„ today would change the epoch of matter radiation 
equality a eq with respect to the massless neutrino case, 
with a scaling a cq oc (1 — fu)^ 1 (assuming ultrarelativis- 
tic neutrinos at the time of equality). This is responsible 
for the main effects on the CMB anisotropy pattern, in 
the range £ < 2.0 eV. Physically, postponing the time of 
equality produces an enhancement of small-scale pertur- 
bations, especially near the first acoustic peak, and in- 
creases slightly the size of the sound horizon at recombi- 
nation. When turning to the growth of structures, there 
is also an additional effect. During matter domination 
and on scales smaller than the free-streaming scale, the 
neutrino perturbations do not contribute to gravitational 
clustering, and neutrinos can be simply omitted from the 
Poisson equation. On the other hand, they do contribute 
to the homogeneous expansion through Friedmann equa- 
tion. Therefore the exact compensation between cluster- 
ing and expansion holding for a pure CDM scenario is 
slightly shifted: the balance is displaced in favor of the 
expansion effect, and the gravitational potential decays 
slowly, while the matter perturbation does not grow as 
fast as the scale factor. Combining the continuity, the 
Euler, the Poisson and the Friedmann equations one gets 



the evolution law for the perturbation in the matter den- 
sity field S m at small-scales 
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which at the first order in the small parameter f v has a 
growing mode solution of the kind S m oc a 1- *^". This is 
valid when p v is dominated by the most massive, non- 
relativistic state(s), and /„ — > constant. Physically, the 
combined effect of the shift in the time of equality and 
of the reduced CDM fluctuation growth during matter 
domination produces an attenuation of perturbations for 
modes k > k nr , where fc nr is the minimum of the comov- 
ing free-streaming wavenumber attained when neutrinos 
turn non-relativistic, and given by [l[ 
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An instantaneous decay of the massive neutrinos at an 
epoch in the matter era can be thought as replacing 
the neutrino fluid with one having the same energy con- 
tent at rjd, but whose energy density scales from that 
moment on as a -4 , since the daughter particles are rel- 
ativistic. Let us estimate how large a value of r/d, or 
equivalently of the proper time td{= r if the neutrino is 
non-relativistic), can be probed cosmologically. Quickly 
after the neutrino decay one has formally f v — > 0, pro- 
vided that td -C to ~ Hq 1 ; from that moment on, the 
cosmological effects of the decaying neutrino scenario are 
analogous to the ones of a massless neutrino universe. 
The condition td *C to is required by the fact that when 
td —> to, the radiation content of the relativistic daugh- 
ters of the massive neutrino has no time to decline to 
zero with respect to the matter density. This condition 
is necessary to change appreciably the energy budget of 
the universe, thus affecting the predicted growth of the 
structures and the time of equality with respect to a mas- 
sive neutrino scenario. Clearly, for a given sensitivity to 
the effect of neutrino masses there is a maximum value 
t™ which would result in a detectable change of cosmo- 
logical observables. A precise estimate of this parameter 
would imply a detailed forecast analysis, which goes be- 
yond the purpose of this paper. Yet, a simple argument 
shows that, relying on the existing forecasts, a conser- 
vative lower limit is £™ ax > rj nr , where i nr is the epoch 
at which the heavier neutrinos become non-relativistic, 
whose redshift is defined by m = 3T„,o(l + z nT ), T v ,o be- 
ing the present temperature of the neutrino gas. Indeed, 
when the decay epoch satisfies td J$ t m , the energy con- 
tent of the products is the same of a relativistic neutrino 
fluid, and it redshifts the same way. So, all physical ef- 
fects of this scenario are basically the same of the case 
where neutrinos are massless 1 . In Fig. [IJ from top to 
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bottom as seen from the left side of the plot, we show 
f v {z) for the following cases: (i) a massive neutrino cos- 
mology, where we assume an IH neutrino mass pattern 
and the lightest neutrino is massless; (ii) as in (i), but for 
NH; (iii) a decaying neutrino cosmology, where massive 
neutrinos have IH; (iv) as in (iii), but for NH; (v) a mass- 
less neutrino cosmology. For the decaying cases, we as- 
sume that all massive neutrinos decay at td = £ nr , where 
i nr is the time of non-relativistic transition of the heav- 
iest neutrino state (m ~ 0.05 eV) . The neutrino mass 
and mixing parameters are from [4|, the cosmic neutrino 
distributions are from jl8j], and for simplicity we have 
assumed a matter-dominated cosmology with the matter 
density parameter flm = 0.24 and the reduced Hubble 
constant h = 0.73 [191 ]. 

Clearly, the cases (iii), (iv), and (v) are very similar 
(exactly degenerate if td <C t ni ) and, as long as t^ < t nT , 
if the massless neutrino case can be disproved, the decay- 
ing neutrino bound immediately follows. The improve- 
ment in the bound on the neutrino lifetime is tremen- 
dous. In particular, neutrinos turn non-relativistic at 
z nl - ~ m/3T U fi ~ 100 mso, i-e. when the universe has 
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of its present age, and 
of the maximum attain- 



able limit reported in Eq. ([Vj, 
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cosmology will be sensitive to E < 0.1 eV, the hierarchy 
information may be available independently. 
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FIG. 1: The function f v {z) for the relevant cosmological cases 
considered in the text. 

To appreciate how strong the bound of Eq. ([5]) would 
be, let us consider a model of a "secret" neutrino interac- 
tion with a (quasi-) massless majoron field <f> of the kind 
C = g DiVj4>+ h.c, i,j labeling different mass eigenstates. 
The total decay rate for a hierarchical neutrino mass pat- 
tern and summing over neutrino and antineutrino final 



state channels is |l4 Il6| 



Obviously, the previous argument does not exclude that 
an accurate forecast analysis may reveal a sensitivity 
to a somewhat larger i™ ax . Note also that we do not 
require that cosmological data need to distinguish be- 
tween NH and IH: if future observations will suggest e.g. 
£ = 0.08 eV with a la error of 0.02 eV, the two neu- 
trino mass patterns would be both consistent within 1 a 
with the best fit, yet a complete decay of neutrinos into 
relativistic particles with lifetime lower than the value 
reported in Eq. ([5]) could be excluded at 4 a. Of course, 
for a given cosmological sensitivity, the significance of the 
above bound increases if the inverted hierarchy is realized 
in nature: in that case £ ~ 0.1 eV holds, and the cosmo- 
logical effects of neutrino masses are larger. Note that 
accelerator neutrino experimen ts l20l| , magnetized detec- 
tors of atmospheric neutrinos [2l|, direct mass searches 
[22] ] , and the serendipitous observation of neutrinos from 
a galactic supernova [23| may all be used to determine 
the mass hierarchy. It is thus possible that by the time 
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are basically insensitive to the energy spectrum of the neutrino 
fluid: They are only sensitive to its overall energy density and 
equation of state. In the case at hand, one or more of the daugh- 
ter particles may have a finite but much smaller mass than the 
parent one, but this does not change our conclusions, at least at 
our level of approximation. 
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This holds in the neutrino rest frame, but in our case this 
is also the Lab decay width, give or take a factor 0(1), 
since the neutrino is just turning non-relativistic. The 
constraint of Eq. ([5]) leads to the stringent bound 



9 <4xl0- 14 m^ 4 . 
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This has to be compared with traditional bounds found 
in the literature in the range g < 10~ 4 -f- 10~ 5 (see e.g. 
(24]]), which is also a typical value invoked in the "neutri- 
noless universe" scenario of Ref . [25[ . Even the extremely 
stringent bound reported in [lt| is more than two orders 
of magnitude weaker. 

Note that the tiny couplings which may induce the 
decay are not sufficient to thermalize extra degrees of 
freedom in the early universe. So, this model does not 
predict departure from the standard expectation for the 
effective number of neutrinos N e ff [l8T | , which can be con- 
sistently fixed in deriving the bound. Yet, if additional 
exotic physics is present, a change (typically an increase) 
of N e g is possible. The effect of a finite £ can be par- 
tially compensated by an increase in N g, which wors- 
ens the sensitivity of cosmological probes to £ (see e.g. 
[2^|). However, the inclusion of these exotic effects in the 
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forecasts can safely follow the standard parameterization 
used in analyses of stable neutrino scenarios. 

Discussion and conclusion Current forecast analy- 
ses suggest that future cosmological surveys may attain 
the sensitivity to detect the effects of a sum of neutrino 
masses as small as ~0.06 eV, the lower limit predicted 
by oscillation data. Provided that the systematics can be 
controlled to that level, we have discussed in this paper 
how such a detection would have profound consequences 
for the particle physics of the neutrino sector, besides 
providing a way to measure the absolute neutrino mass 
scale. In particular, when taking into account the expec- 
tations from the Lab, excluding the £ = case would 
improve by many orders of magnitude the existing limits 
on neutrino lifetime, and as a consequence on neutrino 
secret interactions with (quasi-)massless particles as in 
majoron models. Strictly speaking these bounds apply 
to the heaviest (or the two heaviest, in IH) mass eigen- 
state, but naturaleness and phase-space considerations 
suggest that the lifetime of the lightest state(s) is longer, 
and its coupling with a majoron field weaker, than for 
the heavier one(s). Also, such a bound would be robust 
with respect to the coupling mediating the new inter- 
action (the same may not apply to the considerations of 
[ill, for example). It also applies to any possible invisible 
decay channel, provided that the total mass of the final 
state particles is much smaller than X. In particular, this 
bound applies to 3-v final state decays z/j — ► VjVjVk, as 
well as to decays Vi — > v 3 ■ + <fi in majoron-like models. As 
discussed in [16| , a consequence of such stringent bounds 
is that the decay of high energy neutrinos [28| , a target 
for neutrino telescopes such as IceCube, can not occur. 
Here the conclusion would extend to the diffuse super- 
nova neutrino background, too: a disagreement with as- 
trophysical predictions could not be attribute to neutrino 
decays. We think that the idea developed here provides a 
beautiful example of interplay between particle physics, 
cosmological and astrophysical arguments and motivates 
further the efforts to fully exploit the potential of future 
cosmological surveys. 

Finally, it is worth speculating briefly on the possi- 
bility that, although future observations may attain the 
needed sensitivity, a value of £ consistent with zero is 
favored 2 . This paper suggests that a neutrino lifetime 
t < i nr , as may be due to an extremely tiny coupling of 
the order of<?>4xlCP 14 with a majoron-like particle, 
might provide a possible explanation of a discrepancy 
with oscillation and Lab data, at least if this should arise 
at the S < 0.1 eV level (thus insufficient e.g. to fully ex- 



2 Note that the accuracy needed to detect a finite £ (thus improv- 
ing the bound on the neutrino lifetime) in general differs from the 
accuracy with which £ ~ can be favored: disproving neutrino 
decay may be less challenging than the opposite. 



plain the tension between cosmological mass bounds and 
the claim of detection of neutrinoless double beta decay 
[13] )■ Note that this possibility has some similarity with 
the "neutrinoless universe" scenario of Ref . [IH , since it 
offers a way-out for a possible non-detection of neutrino 
mass in cosmological data, thus re-emphasizing the 
complementarity of cosmological bounds and laboratory 
experiments. However, differently from the latter, it 
would present no departure from the standard cosmology 
as early as the BBN or CMB photon decoupling epoch. 
This avoids completely the constraints based e.g. on N e g 
discussed in [2511 as well as the more stringent arguments 
put forward in [16j ]. 
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